Crystallization of a polymer on a surface 
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We have studied the structure and free energy landscape of a semi-flexible lattice-polymer in the 
presence of a surface. At low temperatures coexistence of two-dimensional integer-folded crystals 
is observed. As the temperature is increased there is a transition from these crystalline config- 
urations to a disordered coil adsorbed onto the surface. The polymer then gradually develops 
three-dimensional character at higher temperatures. We compute the free energy as a function of 
increasing crystallinity and compare with the free energy profiles assumed by the Lauritizen-Hoffman 
surface nucleation theory of polymer crystallization. Our free energy profiles exhibit a 'sawtooth' 
structure associated with the successive formation of chain folds. However, in the early stages of 
crystallization our profiles significantly deviate from those assumed by surface nucleation theory 
because the initial nucleus is not a single stem but two incomplete stems connected by a fold. This 
finding has significant implications for the theoretical description of polymer crystallization. 



I. INTRODUCTION 



On crystallization from solution or the melt simple 
polymers typically form lamellar crystals. As the back- 
bone of the polymer chain is oriented perpendicular to 
the plane of the lamellar, and yet the thickness of the 
crystals are smaller than the length of the chain, a sin- 
gle polymer must traverse the erystal many times folding 
back on itself at each surfaced. Furthermore the thick- 
ness of the crystals have a well-characterized dependence 
on the degree of supercooling. The thickness is always 
slightly larger than the minimum thickness for which a 
lamellar crystal is thermodynamically stable. 

However, although these simple facts were discovered 
about 40 years ago, there is still no consensus on their 
theoretical explanation.^ Moreover, two of the more dom- 
inant apnroaches — Lauritzen-Hoffman surface nucleation 
theoryHTj and the entropic barrier modelETu — appear ir- 
reconcilable. One of the difficulties is that the theories 
often have to make assumptions about the microscopic 
mechanisms of crystallization. As it is difficult to probe 
these processes directly in an experiment, the main test 
for the theories has been the comparison of their pre- 
dictions with macroscopic properties, such as the crystal 
thickness, growth rate and shape. But these tests have 
not proved discriminating enough; the two main theories 
(perhaps with the help of various refinements or 'suitable' 
choices of parameter) have both been able to provide an 
adequate description of many of these properties. 

Therefore, (atomistic) simulation can potentially play 
an important role in this debate by providing insight into 
the microscopic processes involved in polymer crystalliza- 
tion, and thus help in the critical assessment and refine- 
ment of the current theories and perhaps in the devel- 
opment of new th eori es. However, there have been few 
simulation studiest3liil which have attempted to investi- 
gate the basic process in solution crystallization, the ad- 
sorption and crystallization of a polymer on the growth 
surface. (There are many studies which just consider 



adsorption.lI3~fca) Although this is a non-equilibrium pro- 
cess, as a first step it is useful to understand the equi- 
librium behaviour of a single polymer in the presence of 
a surface. In this paper, we first study the basic ther- 
modynamic properties of a simple lattice model of such 
a system (Section HI). This will also add to the increas- 
ing knowledge of the rich phase behaviour of single ho- 
mopolymer chains; in particular, we draw out the simi- 
larities and differences from the behaviour of an isolated 
semi-flexible .polymer, a case that has received much more 
attentionEj~ED Secondly, we examine in detail the free en- 
ergy profile for the crystallization pathway suggested by 
surface nucleation theory (Section |lV|). 



II. METHODS 

A. Polymer Model 

In our model the polymer is represented by an A^-unit 
self-avoiding walk on a simple cubic lattice. There is an 
attractive energy, -e, between non-bonded polymer units 
on adjacent lattice sites and between polymer units and 
the surface, and an energetic penalty, e g , for kinks in the 
chain. The total energy is given by 



'{jlpp -\- Tlp S ^€ -\- TLgCg 



(1) 



where n pp is the number of polymer-polymer contacts, 
n ps is the number of polymer-surface contacts and n g 
is the number of kinks or 'gauche bonds' in the chain, e 
can be considered to be an effective interaction represent- 
ing the combined effects of polymer-polymer, polymer- 
solvent and solvent-solvent interactions, and so our model 
is a simplified representation of a semi-flexible polymer 
at the interface between a polymer crystal and solution. 
The behaviour of the polymer is controlled by the ratio 
fcT/e; large values can be considered as either high tem- 
perature or good solvent conditions, and low values as 
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low temperature or bad solvent conditions. The param- 
eter e g defines the stiffness of the chain. The polymer 
chain is fully flexible at e s =0 and becomes stiffer as e g 
increases. In this study we only consider e g > 0. 

This polymer model has been recently used to study 
the effects of stiffness on .the phase behavipuc-of isolated 
homopolymers by theorycJ and simulatior£J'E£l, and also 
in kinetic Monte Carlo simulations of the growth of poly- 
mer crystals EB 

The global potential energy minimum at a particular 
(positive) e g is determined by a balance between max- 
imizing ripp and rip S , and minimizing n g ; it is a folded 
structure that lies flat on the surface. If the polymer is 
able to form a structure that is a 2-dimensional rectangle 
with dimensions a x b (N — ab), where a < b, then 



= N 

= N -a-b+1 
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have the 



The structures that correspond to n g — r 
polymer chain folded back and forth along the longer 
dimension of the rectangle. By minimizing the resulting 
expression for the energy one finds that the lowest energy 
polymer configuration should have 



a e 



(3) 



Therefore, at e g = the ideal shape is a square and for 
positive e g a rectangle extended in one direction, the as- 
pect ratio of which increases as the chain becomes stiffer. 
However, at most sizes and values of e g it is not pos- 
sible to form a rectangle with the optimal aspect ratio. 
Nevertheless, it is easy to find the global minimum just 
by considering the structures which most closely approx- 
imate this optimal shape. 



B. Simulation Techniques 

To simulate, our system we use configurational-bias 
Monte CarlcO including moves in which a mid-section 
of the chain is regrownEJ We also make occasional bond- 
flipping moves which, although they do not change the 
shape of the volume occupied by the polymp,|jshange the 
path of the polymer through that volume. EStiluj These 
moves speed up equilibration in the dense phases. Dur- 
ing the simulation we always constrain the polymer to 
have one unit in contact with the surface. Thermody- 
namic properties, such as the heat capacity, were calcu- 
lated from the energy distciiuitions of each run using the 
multi-histogram methodBjca 

In order to monitor the state of the polymer we de- 
vised two order parameters. The first, Q2D, probes the 
orientational order within the 2-dimensional polymer ad- 
sorbed onto the surface. 
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where x and y are in the plane of the surface and n a 
is the number of bonds in the direction a. Q2D has a 
value of 1 if all the bonds are in the same direction, i.e. 
the polymer has a linear configuration, and a value of 
if the bonds are oriented isotropically in the plane. The 
second order parameter, Q%d, probes the dimensionality 
of the polymer. 
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Q3D has a value of 1 if all the polymer units are in con- 
tact with the surface and a value of if the bonds are 
oriented isotropically in space. 



In section IV we compute the free energy along path- 
ways characterized by an order parameter, N xta \, which 
measures the degree of crystallinity. This Landau free 
energy is simply related to the canonical probability dis- 
tribution for the order parameter: A^(iV xta i) = A — 
kT log p can (N xta \), where A is the Helmholz free energy. 
However, all relevant regions of this distribution are not 



significantly sampled in t 
we use umbrella sampling 



canonical ensemble, and so 
to calculate the free energy 
accurately over the whole range of the order parameter. 
This is achieved by multiplying the Boltzmann factor by 
the exponential of a biasing distribution, W(N xt ai)', i.e. 
the simulation samples configurations with a probability 
proportional to exp(— (3E + W(N xtEl \)). The canonical 
probability distribution is then obtained from the proba- 
bility distribution from the biased run, Pmulti(iYxtal)j by 

Pcan(N xta l) = Pmulti(Nxtal)exp(-VF(iV x tal))- We wish to 

choose W such that Pmuiti(-^Vxtai) is approximately con- 
stant over the whole. r R n E p of N xta \ (the so-called mul- 
ticanonical approachESEZI) . However, this only occurs 
when W(N xta i) w AL(N xta i)/kT and so we have to con- 
struct W iteratively frorn-the results of a number of short 
preliminary simulations^ 



III. THERMODYNAMICS 



In the presentation of our results we concentrate on 
one example, a 200-unit polymer with e g = At. Results 
for other positive values of e g show the same basic be- 
haviour. We also mainly dwell on those aspects of the 
thermodynamics which are relevant to polymer crystal- 
lization or differ mar kfHJyn from the behaviour of the iso- 
lated homopolymerOci 
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FIG. 1. Some representative configurations for a 200- unit polymer in the presence of a surface, (a) Crystalline configuration 
with 5 stems of length 40 units; it is one of the lowest energy configurations for e 9 = 4e. Q2D=0.960 , i?g=135.3. (b) Disordered 
two-dimensional coil on the surface produced in a run at T = 2.5efc _1 . Q 2D =0.235 , Rg=Hl.6. (c) Three-dimensional coil 
produced in a run at T = 5.0efc _1 . Q 3 d=0.010, ii^ =99.1. (d) Configuration from a simulation at T = 2.375efc~ 1 in which half 
the configuration is crystalline and half is disordered. 



At zero temperature the global potential energy min- 
imum has the lowest free energy. For e g = 4e the op- 
timal aspect ratio of a crystal is 9 (Equation (||)). For 
./V=200 the crystal that most closely approximates this 
shape has five stems (a stem is a straight section of the 
polymer) which are 40 units long; it is the global mini- 
mum. The example shown in Figure [I]a has a tight fold at 
the end of each stem, a situation which is often referred 
to as adjacent re-entry, but there are a number of de- 
generate configurations with non- adjacent re-entry. The 
differences in energy between crystalline configurations 
are only small and so at low temperatures more than 
one type of crystalline configuration is observed in the 
simulation; they can be differentiated by their radius of 
gyration. Figure ^b shows the presence of four coexisting 
crystallites which have four, five, six or seven complete 
stems. This preference for crystalline configuration with 
complete stems is energetic in origin — it maximizes the 
number of polymer-polymer contacts for each fold. A 
similar trend is observed in long monodisperse alkanes; 
the crystals have preferred thicknesses which correspond 



to an integer number of complete stems. I 

As the temperature is increased configurations with 
more disorder in the stem length and a shorter average 
stem length are entropically more favoured. This causes 
the peaks in the probability distribution of the radius of 
gyration to broaden and the radius of gyration to de- 
crease in the range T = 1.0-1. 75efc _1 (Figure ||a). 

As with the isolated homopolymers there must come 
a point when disordered configurations become more 
favoured than crystalline configurations. This 'melting' 
transition is signalled by a heat capacity peak (Figure |3j) 
and by a loss of orientational order (Figure [|a). Although 
Q2D drops down to a value of 0.1 at the transition, Qsd 
remains close to one. This shows that on melting the 
polymer adopts a two-dimensional configuration on the 
surface with no orientational order in the plane (Figure 
|l|b). Furthermore, although this order-disorder transi- 
tion is a finite size analogue of a bulk first-order phase 
transition, no bimodality is seen in the probability dis- 
tribution of Q2D (Figure ^b), i.e. there is no free energy 
barrier between states with high and low Qid- Near to 
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the transition temperature, T m , the Qid probability dis- 
tribution is very broad and flat, perhaps because of the 
contribution of partly crystalline, partly disordered con- 
figurations (such as the one depicted in Figure |l|d) with 
intermediate values of Qid- 



(a) 




^ 2 

FIG. 2. (a) Rg as a function of temperature for a 200-unit 
polymer at e g — 4e. For comparison the results for an iso- 
lated polymer (dashed line) have also been included, (b) 
Probability distribution for Rg at T — 1.25efc -1 . The peaks 
are labelled by the numbers of stems in the corresponding 
configurations. 

The order-disorder transition also causes the radius of 
gyration to increase first as the polymer passes to the 
disordered state and then as this coil begins to expand 
with temperature. However, this rise is checked when 
the polymer begins to develop a three-dimensional char- 
acter (Figure [l|c) . The change in dimensionality occurs 
gradually and is signalled by a decrease in Qsd and also 
a decrease in the radius of gyration to a value which is 
comparable to that for the an isolated chain of the same 
size (Figure ||a). The associated loss in polymer-surface 
contacts is the cause of the high temperature shoulder in 
the heat capacity. Such desorbtion transitions are well- 
understoodliia and we do not dwell on its features here. 
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FIG. 3. C v as a function of temperature for a 200-unit 
polymer at e g — 4e. For comparison the results for an iso- 
lated polymer (dashed line) have also been included. 



In Figures |2j and |3| we have included results for an iso- 
lated polymer for comparison. Until the transition to 
a three-dimensional configuration, the Rg curves for the 
two cases have a similar form (except that the values 
are larger for the 2D polymer on the surface) dua-tp the 
similar nature of the transitions in the two casesrS The 
heat capacity peak for the isolated polymer is sharper, 
suggesting that the melting transition has a stronger 
first-order- like character. Similarly, it is easier to find 
an order parameter for the isolated polymer that shows 
bimodality (i.e. the presence of free energy barrier) in 
the transition region. This difference is probably be- 
cause the energy difference between ordered and disor- 
dered states is less in the presence of the surface; the 
energy of the two-dimensional coil is reduced relative to 
the three-dimensional coil due to the interaction energy 
with the surface. 

The basic behaviour described above generally holds 
for any positive e g , as be can seen from the phase diagram 
in Figure |^. In particular the order-disorder transition al- 
ways leads to a two-dimensional disordered polymer ad- 
sorbed onto the surface, which only at higher temper- 
ature becomes three-dimensional. This fact could be of 
particular relevance to polymer crystallization from solu- 
tion. It suggests that a polymer forms a two-dimensional 
configuration on the growth surface before crystallizing, 
rather than crystallizing directly from solution. The ther- 
modynamics of the crystallization process will depend 
significantly on which of these two possible mechanisms 
holds. The direct mechanism is implicitly or explic- 
itly assumed in several theoretical descriptions E Yet, the 
present results appear to support the assumption made 
by Yamamoto in his study of the dynamics of the crys- 
tallization processJij 
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FIG. 4. (a) Q2D and Q30 as a function of temperature 
for a 200-unit polymer at e g = 4e. (b) Probability distri- 
butions for Q2D at three different temperatures around the 
melting temperature, as labelled. 

However, it is not clear whether this aspect of our 
results — crystallization being preceded by adsorption — 
is likely to be representative of the behaviour of real 
polymers; it might reflect some of the simplicities in the 
model. For example, in our model an adsorbed disor- 
dered polymer and a crystallite can have the same energy 
of interaction with the surface, whereas in reality a poly- 
mer chain that fits snugly into a groove formed by chains 
in the crystalline surface will have a lower energy than a 
random configuration on the surface. Such features could 
result in crystallization and adsorption occurring simul- 
taneously. Therefore, it would be useful if more realistic 
simulations were performed to clarify this situation. 

In the phase diagram in Figure |E| we have not at- 
tempted to distinguish regions where the scaling be- 
haviour of the disordered state is that of a coil or of a 
globule. The main differences in behaviour between poly- 
mers with different values of e g are related to the position 
of the coil-globule transitions with respect to the adsorp- 
tion transition and show up most clearly in quantities 
such as the radius of gyration. For example, the maxi- 
mum in Rg for the two-dimensional disordered polymer 



that occurs at e g — 4e (Figure ||) does not occur for all 
e g . However, we do not pursue these issues further here. 
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FIG. 5. Phase diagram of a 200-unit polymer on a sur- 
face. The phase diagram is divided into regions by the val- 
ues of the melting point, T m , and by the transition from 
two- to three-dimensional configurations, T2-3. These two 
transition temperatures are defined by Q2£>(T m ) = 0.5 and 
Q3d(T2_3) = 0.5. No attempt has been made to differen- 
tiate between coils and globules in the disordered region of 
the phase diagram. 



IV. FREE ENERGY FOR A SINGLE 
CRYSTALLIZATION PATHWAY 

Having understood the basic thermodynamics of our 
system we are now in a position to use the model to ex- 
amine in detail the free energy landscape for the micro- 
scopic mechanism of polymer crystallization assumed by 
the Lauritzen-Hoffman surface nucleation theory.Q~0 In 
this theory the key process in determining the thickness 
of the lamellar crystals is considered to be the nucleation 
and growth of a new polymer layer on the crystal growth 
face (the thin edges of the lamellae). In the theory the 
growth rate for new layers of different thicknesses are 
compared. It is argued that crystals with a thickness 
close to that which gives the maximum growth rate will 
dominate the ensemble of possible crystals, and so the 
predicted thickness of the crystal corresponds approxi- 
mately to the thickness for which the rate is a maximum. 

To calculate the growth rates some assumptions have 
to be made about the processes involved and the associ- 
ated thermodynamics. Firstly, the new layer is assumed 
to grow by the addition of a succession of stems along the 
growth face each connected to the previous by a tight 
fold. Secondly, the length of each stem is assumed to 
be the same as the thickness of the lamella. Thirdly, 
to calculate the free energy for this process macroscopic 
thermodynamic properties are used. This allows the free 
energy as a function of the number of stems to be written 
as 
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A L (N stcnl ) = 2blo + 2(7V stcm - l)abo f - ablAF, (6) 

where o is the surface free energy of the growth surfaces, 
Of is the free energy of the fold surfaces, AF is the free 
energy of crystallization per unit volume, I is the thick- 
ness of the lamella, a is the width of a stem and b is the 
depth of a stem. The geometry of the assumed mecha- 
nism is depicted in the inset of Figure [| The first term 
in the above equation corresponds to the creation of the 
two new lateral surfaces on either side of the nucleus. 
This free energy has to be 'paid for' on the laying down 
of the first stem. The second term corresponds to the 
free energy of the folds that have to be created on the 
deposition of subsequent stems. 

The fourth assumption is that at the barrier between 
states with different N stcm all the new surfaces have been 
created and that a fraction \I> of the free energy of crys- 
tallization has been released. These assumptions lead to 
free energy profiles like that depicted in Figure 0. 




FIG. 6. Free energy profile assumed by the Lau- 
ritzen-Hoffman surface nucleation theory of polymer crys- 
tallization. The theory considers the crystal to form by lay- 
ing down of adjacent stems of the polymer. The inset is 
a schematic representation of a configuration where three 
stems have been deposited. The profile is for a temperature 
at which the crystal is the most stable state 

There have been a number of specific criticisms of this 
free energy profile. Firstly, the model assumes that the 
value of $ for deposition of the first and subsequent stems 
is the same when clearly the physical processes for these 
two cases are very different. The main motivation for this 
assumption is simply that it allows an analytical solution 
to the theory. More recently, there has been some work 
which has attempted to find approximate-solutions when 
this simplifying assumption is relaxedoo 

Secondly, for any non-zero value of 'J the pre- 
dicted thickness becomes infinite at sufficiently large 
supercoolings — the so-called ^/-catastrophe. This oc- 
curs when there is no longer any barrier for the de- 
position of the first stem; i.e. for temperatures where 
$ > 2o/(aAF). To avoid this unwanted behaviour, the 



theory has been supplemented with a justification for a 
zero value of ^. The argument is that the chain first 
forms a weakly physisorbed aligned state which traverses 
the growth face (such a configuration has a high free 
energy because it has lost the entropy associated with 
orientational disorder) before it crystallizes to form the 
first crystalline stemjj rather than the free energy of crys- 
tallization being released as successive units of the first 
stem are deposited onto the surface. There seems little 
evidence to favour the first scenario except that it avoids 
a (^-catastrophe appearing in the theory. Furthermore, 
the second mechanism would be expected to involve a 
lower free energy pathway, and so, if feasible, would be 
preferred. 

If simulations are to help us make some critical as- 
sessment of the free energy profile used in the surface 
nucleation theory we need to devise an order parameter 
which can act as a reaction coordinate along a pathway 
to a specific crystal. We use N xta \ which we define as the 
number of units in the largest fragment of the polymer 
that has part of the structure of our target crystal. This 
target crystal has adjacent re-entry of the stems which 
are are all of a specified length, As our reaction coor- 
dinate is more fine-grained than that used in the surface 
nucleation theory it allows us to find the location of any 
free energy transition states on the path. 

To compare with the theoretical free energy profile the 
pathway should go from the disordered state of the poly- 
mer to the target crystal. However, the difficulty with 
the order parameter, iV x t a i, is that it is not able to dis- 
tinguish between configurations which are disordered and 
those that have a crystalline structure different from the 
target crystal. Therefore, we have to introduce some con- 
straints to ensure that the N — N xts d units which are not 
part of the largest fragment of the target crystal adopt 
a disordered configuration. We use the following order 
parameters to achieve this: Q^di Qj? st an d Qxr- Q^d 
is simply the same as Q2D except that only bonds in 
the non-crystalline part of the polymer are taken into ac- 
count. Q™ st = n T ,? st /2N Tes t, where n T n st is the number of 
instances when a bond in the non-crystalline part of the 
polymer is adjacent and parallel to another bond on the 
surface, and iV res t = N - A'xtal- Qxr = n xr /N xts , h where 
n xr is the number of contacts between the crystalline part 
of the chain and the rest of the polymer. 

The constraints we use are: Q r 2 c ff < 0.3, Q^ st < 0.3 
and Q xr < 0.15. Low values of Q^n an d Q\\ st are appro- 
priate for a disordered chain. (The upper bounds have 
been found to be reasonable from monitoring these order 
parameters for the purely disordered state.) The con- 
straint on Q xr prevents the nucleus being stabilized by a 
more locally ordered region surrounding it. However, the 
exact values for these constraints are slightly arbitrary. 
The free energy profiles that we obtain do reflect these 
choices to a certain extent, but only in the numerical 
details and not in their basic structure. 
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FIG. 7. Free energy for the formation of a target crystal 
with equal stem lengths and adjacent re-entry. In (a) and 
(c) the target crystal has stems 40 units long; the four curves 
correspond to different temperature, as labelled, and we have 
set the zero of free energy to that for the complete crystal. In 
(b) the three curves correspond to target crystals with differ- 
ent stem length, as labelled, all at T = 2.05efc~ 1 ; we have set 
the zero of free energy to that for the complete crystal with 
stem length forty, and the free energy differences between 
the different complete crystals were calculated by complete 
enumeration of possible configurations. The free energy pro- 
files in (a) and (b) are from simulation and those in (c) have 
been calculated using equation ((?]). 



In Figure we show some free energy profiles that we 
obtain for the formation of a target crystal with stems 
that are forty units long at a variety of temperatures. 
It is immediately apparent that the profiles have a saw- 
tooth structure similar to the surface nucleation profiles. 
The minima occur at 21, 31 and Al. Presumably, for longer 
polymers this structure would be repeated at higher mul- 
tiples of I. (Indeed when we use I — 25 minima up to 71 
are observed (Figure 0b).) After the minimum the free 
energy rises sharply with the incorporation of the next 
two polymer units into the crystal reaching a maximum 
at nl + 2 or nl + 3. These steps in the free energy have 
a clear interpretation. Configurations with a complete 
number of stems are favoured because such a configura- 
tion maximizes the number of polymer-polymer contacts 
relative to the number of folds in the crystal. To in- 
crease the size of the crystal fragment a new fold has 
to be formed with an accompanying rise in the free en- 
ergy. After this there is a monotonic decrease in the free 
energy as the new stem grows until the next fold must 
form. Translating this into the language used in the sur- 
face nucleation approach, ^(A^tem — ► Astcm + 1) ~ for 
Astcm > 2 and the transition state occurs virtually im- 
mediately after the completion of the previous stem. 

However, at AT xta j < 21 there are major differences be- 
tween the simulation and the surface nucleation free en- 
ergy profiles. Firstly, there is no sign in the simulation 
results of any feature at iV x tai = I (or at any other value 
of Axtai) due to the formation of the first fold. There- 
fore, the mechanism of its formation must be significantly 
different from that for subsequent folds; it does not oc- 
cur abruptly at a specific value of iVxtal- In the surface 
nucleation approach it is assumed that new stems are 
laid down one at a time. It does not allow for the pos- 
sibility that two new stems could be formed simultane- 
ously. There is no such restriction in our simulations. 
We find that the initial nucleus is not a single stem but 
two incomplete stems of approximately equal length con- 
nected by a fold. As A^tai increases the two stems grow 
in length simultaneously. The reason for this behaviour 
is energetic. It is simple to show that a two-stem nucleus 
can have a lower energy than a single-stem nucleus when 
Axtai > 4e ff /e + 2. Confirmation of this behaviour can 
be found from Figure ||c. The number of gauche bonds 
in the crystal has risen to two by iV x tai ~ 25. The pos- 
sibility of such a twa-stem nucleus has previously been 
suggested by Point O 

We should note that it is possible that the size of the 
surface may have a small effect on the competition be- 
tween the two possible nuclei. On our 'infinite' surface 
the position of the nucleus is irrelevant. However, on 
the growth surface of a lamellar crystal it is important 
that the fold of the two-stem nucleus is close to (or at) 
the edge of the lamella, whereas the initial position of a 
single-stem nucleus is unimportant. Having said this a 
similar preference for two-stem nuclei has been observed 
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in kinetic Monte Carlo simulations of polymer crystal- 
lization on the growth surface of a lamellar crystalES 

Secondly, the behaviour at small iV x tai is different. For 
example at T = 1.7e/c _1 the free energy is downhill until 
the barrier at -/V xta i = 21. For the higher temperature 
free energy profiles there is a free energy minimum corre- 
sponding to the disordered state but it occurs at a finite 
value of -/V x tai, which decreases with increasing tempera- 
ture. At T = 2.05, 2.35 and 2.6ek~ 1 the minima are at 
Nxt&i = 33, 23 and 17, respectively. Furthermore the free 
energy rises very steeply as iV x t a i decreases towards zero. 
These effects occur because crystallization is occurring 
from a disordered state which is adsorbed onto the sur- 
face. This disordered state has structural motifs — single 
straight sections and short folds — which correspond to 
small fragments of the target crystal. For example the 
two-dimensional disordered conformation shown in Fig- 
ure |l|b has -/V x tai = 26. As the temperature decreases the 
persistence length of the polymer increases and so the 
value of JVxtai associated with the coil increases. 

When there is a free energy minimum corresponding 
to the coil, for values of -/V x t a i beyond the minimum the 
free energy rises linearly until iV xta i = 21. Therefore, in 
this temperature range the barrier to forming a crystal 
nucleus with two complete stems does increase with Z, 
as is clearly shown in Figure 0b. However, this is very 
different from the Z-dependence of the initial free energy 
barrier that is predicted by surface nucleation theory; the 
top of the initial surface nucleation barrier occurs near 
to A^tom = 1 (JVxtai = and is a much steeper function 
oil. 

To gain some further insight into the free energy pro- 
files obtained from the simulation, we attempt to model 
them using a simple calculation in which we only explic- 
itly consider the crystalline part of the polymer config- 
uration and assume the rest behaves like an ideal two- 
dimensional coil. The free energy is then 

A(N xtal ) = A coi i(N Iest ) 

(7) 

where the sum is over all possible configurations for the 
crystalline portion of the chain which are N xt sd units long, 
-E x tai , the energy of the crystal, is given by equation (Q) , 
A co n, the free energy of the ideal two-dimensional coil, is 
given by 

^coii(^Vrcst) = -iV res tfcTlog(l + 2exp(-/3e g )) - N TCSt e, 

(8) 

and -Ej i n is the energy of any gauche bonds that might 
by necessity occur at the junction between the crystalline 
and the coil parts of the chain. (The end of an incom- 
plete stem in the crystalline portion of the chain must 
either occur at an end of the chain or be followed by a 
'gauche' bond.) 
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FIG. 8. Comparison of results calculated using equation 
([?]) for crystallization pathways which allow one (dashed line 
1) and two (solid line 2) incomplete stems to simulation (dot- 
ted line s) results, (a) Free energy for the formation of the 
crystal, (b) energy of the crystal (excluding the contribu- 
tion from interactions between the polymer and the surface) , 
E x tai (c) number of 'gauche' bonds in the crystal, 
function of the number of units in the crystal. The crystal 
has stems forty units long. T — 2.05efc~ 1 

Some free energy profiles obtained from this approach 
are shown in Figure Rc. They have a remarkably simi- 
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lar structure to the simulation results, although there are 
two main discrepancies. Firstly, our expression for A co n 
is an underestimate of the true free energy because it ne- 
glects the energetic contribution from contacts between 
different parts of the coil. Therefore, equation (Q) over- 
estimates the stability of the crystal compared to simu- 
lation results at the same temperature. Secondly, we do 
not take into account the fact that structural patterns 
corresponding to small crystal nuclei naturally occur in 
the disordered state, and so the calculated profiles do not 
have a large rise in free energy at small values of AT xtal . 

The useful feature of this simple calculation is that it 
allows us to understand some of the physical origins of 
the features in the free energy profile more easily. For 
example, A(N x t&i — nl) is always significantly lower than 
A(iV x tal = nl — 1) (n > 2); this sharp dip in the free 
energy is because for iV x t a i = nl neither end of the crys- 
talline portion has to be terminated by a 'gauche' bond, 
i.e. Ejoin = 0. 

More interestingly, it can give greater insight into the 
effects of having a two-stem nucleus. In Figure || we 
compare the results when all possible crystalline config- 
urations are considered in the sum of equation (Q) and 
when only those in which there is a single incomplete 
stem are considered. In the latter case the initial nu- 
cleus must be a single stem, and so the results should 
be much closer to the SN profile. They are; the relevant 
free energy profile in Figure now shows a free energy 
barrier for the formation of the first fold at N xt &i — I and 
there is a steep initial rise in the free energy. However, 
it can be clearly seen that the free energy is much lower 
when the possibility of a two stem nucleus is allowed; 
this is because the energy of a sufficiently-long two stem 
nucleus is lower than a single stem (Figure ||b) . Further- 
more, when assuming a single-stem nucleus the results 
for Extai and r7.* tal (Figure ||c) , as well as the free energy, 
significantly differ from simulation results, whereas the 
calculation and simulation are in good agreement when 
the two-stem nucleus is allowed. (The slightly lower value 
of -Extai for the simulation results is due to contacts be- 
tween the crystalline portion and the rest of the chain. 
These are not taken into account in our simple calcula- 
tion.) 

The calculated profiles have two slopes for AT xta i < 2Z. 
For small values of AT xta i the profile is steeper (and al- 
ways positive) because it is still more favourable to have 
a single-stem nucleus. (This feature has no parallel in 
the simulation profile; instead, as discussed earlier there 
is a rapid rise in free energy as iV xta i becomes smaller.) 
In this range the two calculated profiles in Figure ||a are 
very similar. In the range N xta i = 15 — 20 the profile 
changes slope as the structure of the nucleus changes. 
Assuming an ideal two-stem nucleus the slope of the free 
energy profile for 20 < N xt &i < 21 will be 

dA e 

^_ = -- + fcTlog(l + 2e X p(-/3e g ) (9) 
Therefore, the profile in this region is expected to be flat 



for T = 2.03efc~ 1 (Figure ||a) and below this temperature 
the calculated profile will exhibit a maximum at small 
values of A xta i (e.g. T = 1.7efc _1 Figure 0c). A similar 
temperature-dependence of the slope for A xta i ~ 20 — 80 
is seen in the simulation free energy profiles (Figure ^k) . 



V. CONCLUSIONS 

The thermodynamic properties that we find for 
our semi-flexible polymer in the presence of a sur- 
face have-mapy. similarities to an isolated semi-flexible 
polymer jHj'E3€-3n3 e.g. the coexistence of chain-folded 
crystallites with different aspect ratio and the order- 
disorder transition. (It should be remembered that the 
formation of folded structures for our single polymer sys- 
tem is a thermodynamic effect, whereas lamellar poly- 
mer crystals have folded chains due for kinetic, reasons.) 
One new feature is the desorption transitions in which 
the adsorbed two-dimensional polymer gradually adopts 
a three-dimensional configuration. Also of interest is the 
fact that crystallization is always preceded by adsorption. 
However, it is not yet clear whether this feature is likely 
to be common to many polymers or reflects some of the 
aspects of our current model. 

The free energy profiles that we obtain by umbrella 
sampling for specific crystallization pathways have a saw- 
tooth structure where each rise in free energy corre- 
sponds to the formation of a new fold. We found that 
the apportionment factor, ^(JV^em — * Nstem + 1) » for 
A^tcm > 2 (A^stcm is the number of stems in the crystalline 
configuration). However for N stcm < 2 our free energy 
profiles are significantly different from those of Lauritzen 
and Hoffman's surface nucleation theory. There is no 
feature in our profiles that corresponds to the formation 
of the first fold because the initial nucleus is not a sin- 
gle stem but two incomplete stems connected by a single 
fold. The nucleus increases in size by the simultaneous 
growth of both stems. The reason that a two-stem nu- 
cleus is preferred is simply because beyond a certain size 
it has a lower energy than a single-stem nucleus. Such 
energetic considerations are also likely to hold for more 
realistic polymer models. This finding has serious impli- 
cations for the coherence of the surface nucleation theory 
of polymer crystallization because the barrier to the for- 
mation of the first stem, and its dependence on crystal 
thickness, plays a key role in the theory. 

Our results also hint at another potential problem with 
the surface nucleation theory To obtain the free energy 
profiles of Figure ^ we had to constrain the system to pre- 
vent crystallization occurring by other pathways. In sur- 
face nucleation theory only a small subset (those where 
the stem length is constant) of the multitude of possible 
pathways are taken into account. However, a rigorous 
treatment pnuires that the full web of possible paths 
is explored OC3 Recent kinetic Monte Carlo simulations 
which do just this present a radically different picture of 







the mechanism by-which the thickness of lamellar crys- 
tals is determined H3 
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